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Abstract. Self-similarity in general relativity is briefly reviewed and the differences between self- 
similarity of the first kind (which can be obtained from dimensional considerations and is invariantly 
characterized by the existence of a homothetic vector in perfect fluid spacetimes) and generalized self- 
similarity are discussed. The covariant notion of a kinematic self-similarity in the context of relativistic 
fluid mechanics is defined. It has been argued that kinematic self-similarity is an appropriate generaliza- 
tion of homothety and is the natural relativistic counterpart of self-similarity of the more general second 
(and zeroth) kind. Various mathematical and physical properties of spacetimes admitting a kinematic 
self-similarity are discussed. The governing equations for perfect fluid cosmological models are intro- 
duced and a set of integrability conditions for the existence of a proper kinematic self-similarity in these 
models is derived. Exact solutions of the irrotational perfect fluid Einstein field equations admitting a 
kinematic self-similarity are then sought in a number of special cases, and it is found that; (1) in the 
geodesic case the 3-spaces orthogonal to the fluid velocity vector are necessarily Ricci-fiat and (ii) in 
the further specialisation to dust (i.e., zero pressure) the differential equation governing the expansion 
can be completely integrated and the asymptotic properties of these solutions can be determined, (iii) 
the solutions in the case of zero-expansion consist of a class of shear-free and static models and a class 
of stiff perfect fluid (and non-static) models, and (iv) solutions in which the kinematic self-similar vec- 
tor is parallel to the fluid velocity vector are necessarily Friedmann-Robertson- Walker (FRW) models. 
Solutions in which the kinematic self-similarity is orthogonal to the velocity vector are also considered. 
In addition, the existence of kinematic self-similarities in FRW spacetimes is comprehensively studied. 
It is known that there are a variety of circumstances in general relativity in which self-similar models 
act as asymptotic states of more general models. Finally, the questions of under what conditions are 
models which admit a proper kinematic self-similarity asymptotic to an exact homothetic solution and 
under what conditions are the asymptotic states of cosmological models represented by exact solutions 
of Einstein's field equations which admit a generalized self-similarity are addressed. 



1. Introduction 

Self-similar solutions were originally of interest since the governing equations of a given problem 
simplify and often systems of partial differential equations (PDEs) reduce to ordinary differential 
equations. Indeed, self-similarity in the broadest (Lie) sense refers to an invariance which simply 
allows the reduction of a system of PDEs. Self-similarity refers to the fact that the spatial distribution 
of the characteristics of motion remains similar to itself at all times during the motion and self-similar 
solutions represent solutions of degenerate problems in which all dimensional constant parameters 
entering the initial and boundary conditions vanish or become infinite (Barenblatt and Zeldovich, 
1972). Indeed, such solutions describe the "intermediate-asymptotic" behaviour of solutions in 
the region in which a solution no longer depends on the details of the initial and/or boundary 
conditions. Cases in which the form of the self-similar asymptotes can be obtained from dimensional 
considerations are referred to as self-similar solutions of the first kind (Barenblatt and Zeldovich, 
1972). 

Similarity solutions are of importance within general relativity. For example, a strong explosion in 
a homogeneous background produces fluctuations which may be very complicated initially, but they 
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tend to be described more and more closely by a spherically symmetric similarity solution as time 
evolves (Sedov, 1967), and this applies even if the explosion occurs in an expanding cosmological 
background (Schwartz et al., 1975; Ikeuchi et al., 1983). The evolution of voids is also described by 
similarity solutions at late times (Bcrtschingcr, 1985). In addition, the expansion of the universe 
from the big bang and the collapse of a star to a singularity might (both) exhibit self-similarity in 
some form since it might be expected that the initial conditions 'are forgotten' in some sense. 

In this paper we shall assume that the source of the gravitational field is that of a perfect fluid; 
i.e., the energy- momentum tensor is given by 

Tab = (/U + p)u a Ub + pg a b, (1-1) 

where u a is the normalized fluid 4-velocity and /z and p are, respectively, the density and pressure. 
In natural units c = 8nG = 1, the Einstein field equations (EFEs) of general relativity then read 

G a b = T ab , (1-2) 

where G a b is the Einstein tensor. 

Similarity solutions were first studied in general relativity by Cahill and Taub (1971), who did so 
in the cosmological context and under the assumption of a spherically symmetric distribution of a 
sclf-gravitating perfect fluid. They assumed that the solution was such that the dependent variables 
are essentially functions of a single independent variable constructed as a dimcnsionlcss combination 
of the independent variables and that the model contains no other dimensional constants. Cahill 
and Taub (1971) showed that the existence of a similarity of the first kind in this situation could be 
invariantly formulated in terms of the existence of a homothetic vector. A proper homothetic vector 
(HV) is a vector field £ which satisfies (after a constant rescaling) 

£tg a b = 2g ab , (1.3) 

where g a b is the metric and C denotes Lie differentiation along £. 
It follows from equation (1.3) that 

C^Gab = 0. (1.4) 

When the source of the gravitational field is a perfect fluid, in the case of self-similarity of the first 
kind it follows from dimensional considerations that the physical quantities transform according to 

L iU a = -u a , (1.5) 

and 

C^fj, = —2[i, C^p = —2p. (1.6) 

From these equations it follows that 

C^Tab = 0, (1.7) 

which is therefore consistent with the EFEs (1.2). Indeed, in the case of a perfect fluid it follows that 
equations (1.5) and (1.6) result from equations (1.3) [through eqns. (1.2), (1.4) and (1.7)] so that 
the physical quantities transform appropriately (Cahill and Taub, 1971; Eardley, 1974). Hence in 
this case "geometric" self-similarity and "physical" self-similarity coincide. However, this need not 
be the case (see Coley, 1996, for details). The properties of the matter and those of the geometry 
are related through the EFEs, and in general there will be further constraints arising from the 
compatibility of the EFEs and the conditions of self-similarity ( "integrability" conditions). We shall 
investigate these integrability conditions later. 

A. Self-similarity of the second kind. 

The existence of self-similar solutions of the first type is related to the conservation laws and to 
the invariance of the problem with respect to the group of similarity transformations of quantities 
with independent dimensions, in which case a certain regularity of the limiting process in passing 
from the original non-self-similar regime to the self-similar regime is implicitly assumed. However, 
in general such a passage to this limit need not be regular, whence the expressions for the self- 
similar variables are not determined from dimensional analysis of the problem alone. Solutions 
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are then called self-similar solutions of the second type. Characteristic of these solutions is that 
they contain dimensional constants that are not determined from the conservation laws (but can be 
found by matching the self-similar solutions with the non-self-similar solutions whose asymptotes 
they represent) (Barcnblatt and Zeldovich, 1972). 

Self-similarity in the broadest sense refers to the general situation in which a system is not 
restricted to be strictly invariant under the relevant group action, but merely to be appropriately 
rescaled. The basic condition characterizing a manifold vector field £ as a self-similar generator 
(Carter and Henriksen,1991) is that there exist constants di such that for each independent physical 
field 

Ce*\ = di*\, (1.8) 

where the fields $^ can be scalar (e.g., ji), vectorial (e.g., u a ) or tensorial (e.g., g a b)- In general 
relativity the gravitational field is represented by the metric tensor g ab , and an appropriate definition 
of "geometrical" self-similarity is necessary. In the seminal work by Cahill and Taub (1971) the 
simplest generalization was effected whereby the metric itself satisfies an equation of the form (1.8), 
namely £ is a HV, this evidently corresponding to Zeldovich's similarity of the first kind. 

However, in relativity it is not the energy-momentum tensor itself that must satisfy (1.8), but each 
of the physical fields making up the energy-momentum tensor must separately satisfy an equation 
of the form (1.8). In the case of a fluid characterized by the timelike congruence u a the energy- 
momentum tensor can be uniquely decomposed with respect to u a (Ellis, 1971), and each of these 
uniquely defined components (each of which has a physical interpretation in terms of the energy, 
pressure, heat flow and anisotropic stress as measured by an observer comoving with the fluid) 
must separately satisfy an equation of the form (1.8). In the same way, if the metric can be 
uniquely physically and covariantly decomposed then the homothetic condition can be replaced by 
the conditions that each uniquely defined component must satisfy (1.8), maintaining self-similarity 
For example, in the case of a fluid, the metric can be decomposed uniquely in terms of u a , through 
the projection tensor 

Kb= 9ab + U a U b , (1.9) 

into parts h a b and (minus) u a Ub- The projection tensor represents the projection of the metric into 
the 3-spaces orthogonal to u a (i.e., into the rest frame of the comoving observers), and if u a is 
irrotational these 3-spaces are surface forming, the decomposition is global, and h a b represents the 
intrinsic metric of these 3-spaces. [h a b is the first fundamental form of the hypersurfaces orthogonal 
to u a and can be regarded as the relativistic counterpart of the Newtonian metric tensor, when 
the flow independent u a is defined as the relativistic counterpart of the preferred ( irrotational ) 
Newtonian time covector — t a (Carter and Henriksen, 1989).] 

B. Kinematic self-similarity. 

It is arguments similar to these, and, more importantly, a detailed comparison with self-similarity 
in a continuous Newtonian medium, that has led Carter and Henriksen (1989) to the covariant notion 
of kinematic self-similarity in the context of relativistic fluid mechanics. A kinematic self-similarity 
vector £ satisfies the condition 

C(U a = au a , (1-10) 

where a is a constant (i.e., £ is a continuous kinematic self-similar generator with respect to the flow 
u a ). Furthermore, 

C ( h ab = 2h ab (1.11) 

(£ has been normalized so that the constant in (1.11) has been set to unity). Evidently, in the case 
a = 1 it follows that £ is a HV (Cahill and Taub, 1971), corresponding to self-similarity of the 
first kind (Barrenblatt and Zeldovich, 1972). Carter and Henriksen (1989) then argue that the case 
a 7^ 1 (a = 0) is the natural relativistic counterpart of self-similarity of the more general second 
kind (zeroth kind). 

The parameter a represents the constant relative proportionality factor governing the rates of 
dilation of spatial length scale and amplification of time scale. Evidently, when a^l (i.e., £ is not 
a HV), the relative rescaling of space and time (under £) are not the same (and in the zeroth case 
there is a space dilation without any time amplification). 
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C. Remarks. 

1. All of the fields must satisfy equations of the form (1.8). However, not all of the di need be 
independent. In practice it is convenient to assume independence and then determine the constraints 
on the di (in addition to the constraints on the form of the solutions) that arise from the imposition 
of the equations of state and the EFEs. In the case of a perfect fluid we have that, in addition to the 
conditions of "kinematic" self-similarity and "geometric" self-similarity, as represented by equations 
(1.10) and (1.11), respectively, the further conditions of "physical" self-similarity, given by 

jC^fj, = a^i, £^p~bp, (1-12) 

must be satisfied, where a and b are constants. That is, in order for a perfect fluid spacetimc to 
admit a proper kinematic self-similarity all of the conditions represented by equations (1.10)— (1.12) 
must be satisfied. 

2. Additional constraints may arise from the imposition of an equation of state. The least 
restrictive case arises from the exceptional pressure-free case. Here p = 0, whence 

C^T ab = (2a + a)T abl (1.13) 

and no further restrictions apply. Usually a linear barotropic equation of state of the form 

P=(7-1)M, (1-14) 

is assumed, where the constant 7 obeys 1 < 7 < 2 for ordinary matter. [The case < 7 < 2/3 is 
of interest in studying models that undergo inflation.] If a barotropic equation of state p = p(/j.) is 
assumed, it is well known that if a spacetime admits a non-trivial HV then equation (1.14) necessarily 
results (Cahill and Taub, 1971). If, on the other hand, a spacetime admits a non-trivial kinematic 
self-similarity then (except in the special case of dust) the polytropic equation of state p = Pol^ 1 
follows from equations (1.12) (where 7 is the polytropic index and 7 = b/a). 

We note from equation (1.13) that in the case of dust the total energy-momentum tensor T ab 
itself trivially satisfies an equation of the form (1.8). When p ^ 0, this is only be possible in the 
special case that b = a + 2(a — 1). 

3. The differential geometric properties of HV were studied by Yano (1955). The totality of 
HV on a spacetime form a Lie algebra H n (of dimension n) which contains (if H n is non-trivial) 
an (n — 1) dimensional (sub)algebra of KV, G n -i. If a given spacetime is not conformal to an 
exceptional 'plane wave spacetime', it follows that if the orbits of H n are r-dimcnsional, then the 
orbits of G n -i are (r — l)-dimcnsional (Eardley, 1974), and if a spacetime is not conformally flat, 
the spacetime is conformally related to a spacetime for which the Lie algebra H n is the Lie algebra 
of KV (Defrise-Carter, 1975). 

The totality of kinematic self-similar vector fields on a spacetime also form a Lie algebra (Carter 
and Hcnrikscn, 1991) which wc shall denote here by K. Now, H C K, but K neither contains nor 
is contained within either the conformal algebra C in general or the inheriting algebra / (Castejon- 
Amenedo and Coley, 1992) in particular. 

Let us suppose that £1 and £ 2 € K, so that 

£t 1 h ab = 2h ab , C^Ua^aiUa, 
C^ 2 h ab = 2h ab , C^ 2 u a = a 2 u a , 

where £1 and £2 are linearly independent and «i / 1 and a 2 7^ 1 
defining 

£ = (4i+d£ 2 , 
where c and d are non-zero constants, we find that 

C i h ab = 2(c + d)h ab (1.16) 

and 

£ € u a = (cat + da 2 )u a . (1-17) 



(i.e., neither are HV). Then, by 

(1.15) 
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[We note that equations of the form (1.12) are trivially satisfied for £ given by (1.15)] . Now, if u\ ^ a 2 
(where either ol\ or a 2 may be zero), then we can always choose d = 1 — c and c = (1 — a 2 )/(a 1 — a 2 ) 
so that C^h a b — 2h a b and C^u a = u a ; that is, £ is a HV. If, on the other hand, a\ — a 2 (possibly 
zero), then we can choose c + d = 0, whence C^h a b = and C^u a = so that £ is, in fact, a 
KV. Therefore, we have shown that each non-trivial K n {n > 1) contains an (n — l)-dimcnsional 
(sub)algebra of HV, H n _i, where H n -\ m &y be trivial (i.e., H n _\ need not contain any proper HV). 
An illustration of this result can be found in section 4.C. 

4. In the case of radiation (i.e., p = so that T = 0) the existence of a HV implies the existence 
of a conserved quantity. In general, if we define the current P a — T ab £b, then from energy- momentum 
conversation, P a - a = T ab g ab = T. For radiation, P a = § (4u a [u fc £ b ] +£ a ), and P a , a = 0. In the case 
of kinematic self-similarity it follows from equations (1.10) and (1.11) that P a - a = 3p—a^i, implying 
the existence of a conserved quantity for fluids with an equation of state p = f/i. In particular, in 
the case of kinematic self-similarity of the zeroth kind (i.e., a = 0), there exists a conserved quantity 
in the special case of dust (i.e., p = 0). 

5. Finally, there is another case of potential interest in which the ratio a/1 (where a is defined in 

(1.10) and we recall that the constant in (1.11) has been normalized to unity) approaches infinity, 
and we could refer to this case as kinematic self-similarity of 'infinite' kind. This case could be 
covariantly defined by equation (1.10) (in which a could be normalized to unity) and equation 

(1.11) with zero right-hand side (i.e., C^h a b = 0), and £ consequently represents a generalized "rigid 
motion" . This case will be investigated elsewhere. 

In section 2 we shall describe the cosmological models under investigation and introduce their 
governing equations. We shall then obtain the set of integrability conditions for the existence of a 
kinematic self-similarity in the models. Finally, we shall deduce the equations in the relevant case 
of zero vorticity. In sections 3 and 4 we shall study the models in a number of special cases; namely, 
in the case of zero acceleration (3. A) and the further subcase of zero pressure (3.B), the case of 
zero expansion (which includes the special subcase of static models) (3.C), the cases in which the 
kinematic self-similarity is cither parallel to or orthogonal to the fluid velocity vector (4. A and 4.B, 
respectively) and finally we shall investigate the existence of kinematic self-similar vectors in FRW 
spacetimes. In section 5 we shall summarize the main results and we shall discuss further avenues 
of research. 

2. Analysis 

A. The governing equations. 

The covariant derivative of u a can be decomposed according to (Ellis, 1971) t 

U a; b = (Tab + ^Ohab + U a b - U a Ub, (2.1) 

where B a b = h( a c h b f u C :dj = g ab a b is the expansion, a a b = B a b — \®h a b is the shear tensor and 
cr 2 = \o a b<J ah ', uj a b = h[ a c hbf u c;d is the vorticity tensor and uj 2 = ^uj a buj ab , and u a = u a -bu b is the 
acceleration.''' Using these definitions the governing equations can be written down (Ellis, 1971). 
The conservation laws, in the case of a perfect fluid, become 

fi+(p + p)0 = O, (2.2) 
(n+p)u a = -PAh b a , (2.3) 

where fi = n^ a u a (for example). Equations (4.12), (4.15), (4.17), (4.18), (4.16) and (4.19) in Ellis 



tWe shall follow the notation and conventions in Ellis (1971); in particular, Roman indices range from to 3 and 
Greek indices from 1 to 3. 



(1971), become, respectively, 



6 + \o 2 - u a ia + 2(<r 2 - u 2 ) + + 3p) = 0, (2.4) 

h a b {^P^j = &° b (u b + pcu b ^j = a a b J> + \rf bcd u b u c , d , (2.5) 

h e b (u bc ;c - a bc , c + + (co e 6 + a e fc )7i b = 0, (2.6) 

w a , a = 2uj b ii b , (2.7) 

fib* h b 9 frf g - ft, a / /l b 9 U(f-g) - U a U b + lO a U b + (JafCF 1 b 

+ ^9<i ab + h ab (^-^uj 2 - 2 -o 2 + + £ ab = 0, (2.8) 

H ad = 2u {a w d) - hj h d °(uj ig h;c + a {g b;c )v s )fbcU f , (2.9) 

where E ab and H ab are, respectively, the electric and magnetic parts of the Weyl tensor.^ 

The remaining equations are the Bianchi identities (cf. Ellis, 1971, eqns. (4.21a -d)) and the field 
equations (cf. Ellis, 1971, eqns (4.23) - (4.26)). In particular, in the case of zero vorticity (w = 0) 
the Gauss-Codacci equations become 

3 Rab = {-0(7 f g - &fg + U(f;g))hJ h b 9 

+ u a ii b + ^h ab (-^9 2 + 2a 2 + 2n-ii c , c ^j , (2.10) 

and 

3 R = -?-6 2 + 2a 2 + 2^, (2.11) 
o 

where 3 R ab is the Ricci tensor of the 3-spaces orthogonal to u a and 3 R is the corresponding Ricci 
scalar. 

B. Integrability Conditions. 

From equations (1.10) and (1.11), the existence of a kinematic self-similar vector implies that 

^9ab = 2g ab + 2(1 - a)u a u b , (2.12) 
where a ^ 1 is assumed hereafter. From Yano (1955) we have that 

CtR ab = (£ 5 r^) ;c - (C^T c ac ), b (2.13) 

where 

W c = \g at [{L^ bt ), c + (jC i9ct ) ; b - [C^bch]- (2-14) 

Hence we have that 

C^R ab = (1 - a)g ct [(u b u t )-ac + {u a u t )- bc - {u b u a ) ;t c\- (2-15) 
Therefore, using (2.1), after a long calculation we obtain 
1 



(1-a) 



C^R ab = 2a ab + 29a ab + 2a bc uj c a + 2a ac uj c b 
2 

+ 4uj ac uj c b + -9g ab [9 + 9}+ u a [2ui cb ' c - 2a bc ii c - 2u bc u c 



+ u b [2uj C a ' c - 2<7 ac u c - 2ui ac u c ] + u a u b 



^9{9 + 9)-2u c - c 



(2.16) 



Thus, decomposing equation (2.16) using u a , we obtain 
1 



— -£^R a b 



1 



L(l-a) 
1 



C^Rab 



(1-a) 
1 — a 



C-iRab 



u a u b = -8uj 2 - 2u c ;c , 
fr afc = 2(0 + 2 -4a; 2 ), 
u a /i^ = 2w da u a + 2w dc ;c - 4w 2 u d , 



1 



^ - -h ef h' 



ah 



2[& e f - UfUebU - U e (7fbU 



+ 6a ef + Gf c uj° e + a ec u)°f + 2u e c u) cf + -h ef uj 2 }. 



In the case of a perfect fluid, from the EFEs we have that 

Rab = ^(n + 3p)u a u b + ^(fi-p)h ab , 
whence, using equations (1.10) and (1.11), and (1.12), we obtain 



C^Rab = -^{{a + 2a) fj, + 3(6 + 2a)p}u a u b 
+ I{( a + 2)/i- (b + 2)p}h ab . 



Decomposing equation (2.22) yields 

1 



(1-a) 
1 



[C^Rab] U a U b = 



a + 2a 



2(1- 


a) 


3(a + 


2) 


2(1- 


a) 



/i + 
/i + 



3 [6 + 2a] 

-3(6+2)' 
2(1 -a) _ 



(2.17) 
(2.18) 
(2.19) 

(2.20) 
(2.21) 



(2.22) 



(2.23) 



(1-a) 

[C^Rab] u a h b d = 0, 

[C^Rab] (h a e h b f - \h ef h ab ^j = 0. 

Hence, the integrability conditions for the existence of a proper kinematic self-similarity become 



(2.24) 
(2.25) 
(2.26) 



- 8lu 2 - 2u c ;c = 
2(6 + 2 -4uj 2 ) = 



a + 2a 
2(1 - a) 
3(a + 2) 



2(1 -a) 
2u; da u a + 2uj dc ;c - 4cj 2 u d = 0, 



3[6 + 2a] 
2(1 -a) 
3(6 + 2) 



2(1 -a) 



(2.27) 

(2.28) 
(2.29) 



and 



& e f-Ufa eb u -u e <T fb u + 9a ef + a fc uj c e + 

<J ec UJ C f + 2uJ e C Ld cf + ^h ef 0J 2 =0. 



(2.30) 



Equations (2.27)-(2.30) must be satisfied in addition to equations (2.2)-(2.9). 
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C. Vorticity-free case. 

In the Introduction we discussed the analogy with self-similarity in the Newtonian case where the 
irrotational case was singled out to be of special importance. Setting ui = in the above equations 
we obtain: 



fi+ (p+p)9 = 

(fj, +p)u a 



-p,bK, 

<+^9 2 -u a . a + 2a 2 + ^(p + 3p) = Q, 



and 



-ct oc ;c + -^ ) h% + a%u b 



0. 



The integrability conditions (2.27)-(2.30) become 

- 2u c ;c = 

2(6 + 6 2 ) = 



" a + 2a ' 




"3(6- 


-2a)' 


[2(1 - a) _ 


[2(1 


— a) 



p, 



"3(a + 2)" 


PL + 


"-3(6- 


1-2)1 


[2(1 - a) _ 


L 2(1- 


a)\ 



p, 



(2.31) 
(2.32) 

(2.33) 
(2.34) 

(2.35) 
(2.36) 



and 



& e f - Uf(T e bU b - U e (TfbU b + 9a ef = 0. 



(2.37) 

Finally equations (2.10) and (2.11) must be satisfied. The remaining non-trivial equations are 
equations (2.8) and (2.9) which serve to define E a b and H a b, respectively We note that from 
equations (2.8), (2.10), (2.11) and (2.37) wc have that 



1 



'Rab - 7;h ab 3 R = E ab + cr a /(T J f h — -8a ab — -9 2 h a b 



Now, contracting equation (2.37), we obtain 



Using (2.35), equation (2.33) becomes 

3 4(1 — a) 



3(6 + 2) 



p = 0. 



4(1 -a) 1 

Differentiating this expression, and using equations (2.39) and (2.40), yields 



9 + 



+ 



+ 



(a + 2) 



[A + 20/i] + 



3(6 + 2) 



[p + 29p] = 0. 



(2.38) 
(2.39) 
(2.40) 

(2.41) 



3"" ' 3" ' 4(1 -a) ^ ^' 4(1 - a) 
Finally, using equations (2.31) and (2.36), we obtain (9 ^ 0) 

* + * 2 = 4 {r^- rt - (2 - 42) 

Hence from (2.36) we have that 

(a - b)p = 0. (2.43) 

We shall study the special case of dust (p = 0) in the next subsection. If p ^ 0, then necessarily 
a = b. 

In either case, equations (2.40) and (2.42) yield 



-^9 2 + 2a 2 



From equation (2.11) we then have that 



, (Q + 2) 
(1-cO 

(a + 2a) 



H = 0. 



(2.44) 



(2.45) 



(I -a) 

From the above equations we note a number of interesting special cases arising. When p ^ 0, we 
see that if p = [i (stiff matter) ora = —2, equation (2.42) reduces to a simple first order DE for 9. 
In addition, if a = —2, then from (2.44) we see that 9 2 = 3a 2 . Finally, if a = —2a, then from (2.45) 
we see that the Ricci scalar curvature of the 3-spaces orthogonal to u a vanishes. 
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3. Special Cases 

A. Geodesic case. 

Let us first consider the case in which the acceleration u a is zero. First, we note that when u a is 
irrotational and geodesic there exists preferred coordinates in which (Coley and McManus, 1994) 

ds 2 = -dt 2 + H a/3 {t,x"')dx a dx 13 , (3.1) 

and 

u a = -6° a (3.2) 

(and in these coordinates we now have, for example, 9 = 8, t )- We also note that if the shear 
is zero, a = 0, then the spacetime is necessarily FRW (Ellis, 1971; Coley and McManus, 1994). 
Henceforward we shall assume that a ^ (FRW spacetimes will be considered later). 
Immediately, from (2.32) we have that 

P, b h b a = (3.3) 
(or p tCl = in the preferred coordinates), and since 

(a-b)p = 0, (3.4) 

equation (2.35) yields (/z + 3p ^ 0) 

a = -2a, (3.5) 
whence (in both cases a = b and p = 0) it follows from equation (2.42) that 

+ 9 2 -^+^ = 0, (3.6) 
and consequently from (2.33) or (2.44) we obtain 

-^6> 2 + 2cr 2 + 2/i = 0. (3.7) 
Therefore, from (2.11) or (2.45), we find that 

a R = 0, (3.8) 

and since, from equation (2.37), we have that 

& ef + 6a ef = 0, (3.9) 

and hence 

0- + 6>o- = 0, (3.10) 
from equation (2.10) we consequently have that 

3 R ab = 0; (3.11) 

i.e., the 3-spaces orthogonal to u a are, in fact, Ricci flat. 

Spacetimes in which the 3-spaces orthogonal to u are Ricci flat have been studied by a number of 
authors (e.g., Collins and Szafron, 1979; Stcphani and Wolf, 1986; see also Kramer et al., 1980). In 
particular, examples in which u is irrotational and geodesic include a subclass of orthogonal spatially 
homogeneous models (i.e., the Bianchi I spacetimes; see Kramer et al., 1980) and a subclass of the 
Szekeres cosmological models (Collins and Szafron, 1979; Goode and Wainwright, 1982). We shall 
return to this and attempt to exploit previous work elsewhere. 
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The remaining non-trivial equations are [(2.31) and (2.34)] 

A + {jJL + p)e = o, (3.12) 

and 



Hr b % + V> e 6 = 0, (3.13) 



which implies that 



, 2 • 
2cr 2 + -6 



u e , (3.14) 



and equations (2.8) and (2.9) yield simplified expressions for E a b and H^. Finally, differentiating 
(3.6), and using equations (3.6) and (3.12), we obtain 

9 + 390 + 9 3 = -|(p + 2p6»). (3.15) 

In the preferred coordinates (3.1), equations (1.10) and (3.2) yield 

C, a = a6Z; £° = at + c, (3.16) 

and equations (1.12) and (3.3) yield 

£<PW) = g®=, (3-17) 
at at + c 

the solution of which can be written as (a ^ 0^) 

P = Pot- 2 , (3.18) 

since from (3.5) either b = —2a or p = 0, and we have set c = so that p — > oo as t — > + . In 
addition 1 *, equation (1.11) yields 

r,* = o ; r = r(^ 7 ), (3.19) 

and from equations (1.12) and (3.12) we obtain 

^ jQ £ a = -lap, + a( Po + pt 2 )9t~ x . (3.20) 
Next, the differential equation (3.15) becomes 

Q.tt + 309, t + 3 = 3 Po t- 3 {l - 9t). (3.21) 
Defining the new time, t, and the new variable, ip, by 

^=0, ^ = 9\ (3.22) 

equation (3.21) becomes 

ip" + ty' + 2?A = 3 Po t- 3 (ip- 1/2 - t), (3.23) 

where ip' = f£ and 

t = J tp-^ 2 dT. (3.24) 



tThe special case a = will be dealt with separately at the end of this subsection. 
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In principle we can solve equations (3.23) and (3.24) for ip — ip{r, x 1 ) and r = r{t, x 1 ) in terms 
of arbitrary functions of x 1 . We can then integrate equation (3.10) to obtain 

a 2 = £ 2 (xT)e- 2T . (3.25) 

Equation (3.9) becomes 

2 1 

a 'afi - -Q a aiVp 7 + gCr Q/3 = 0. (3.26) 

In addition, equation (3.7) yields 

/i = -£ 2 e- 2T + ^V, (3.27) 

where (eqn. (3.12)) 

n' + H=-p t- 2 . (3.28) 

The remaining non-trivial equations serve to determine the metric functions H a p [e.g., cqns. 
(1.11) and (3.11)], or to constrain (through differential relationships) the various arbitrary functions 
of x"< [e.g., eqns. (3.14), (3.20), (3.27) and (3.28)]. The special case p = will be dealt with in the 
next subsection. If p ^ and there exists an equation of state of the form p = p(p), then from 
equations (1.12) and (3.4) we obtain p — (7 — and hence 

fi = ^i t- 2 . (3.29) 

From (3.12) we then obtain (fi + p a ^= 0) 

°- 2M ° r\ (3.30) 



Vo+Po 

whence equation (3.6) yields 

2 

Mo =Po or p = -^=^o /2 -Mo- (3-31) 

In the latter case 8 2 = 3/i whence from (3.7) a — and the spacetime is necessarily FRW, and in 
the former case we have the special case of a stiff fluid with 8 — t^ 1 and a 2 = ii~ 2 (l — 3/i ) [and, of 
course, 3 i? a {, = ; eqn. (3.20) is satisfied identically and a a b satisfies eqns. (3.9) and (3.25), where 
e~ T = t -1 , and a a c ;c = 0]. Finally, we note the special non-trivial solution of Benoit and Coley 
(1996) in the particular case of spherical symmetry. 

Special case a = 0. 

In the special case of kinematic self-similarity of the zeroth kind (i.e., a = 0), the analysis up to 
equation (3.19) is similar resulting in 

a = b = Q, (3.32) 

and 

p = p \ p constant , (3.33) 

where 

r = (c,r(* 7 )). (3.34) 

The differential equation (3.23), obtained from (3.15), becomes 

ip" + + lip = -Q Po . (3.35) 

This equation can be integrated to obtain 

9 2 = V = c(x a )e- 2T + al(x a )e- T - 3 Po , (3.36) 

and the analysis then essentially follows the same steps as in the zero-pressure case (see the details 
in the next subsection). Indeed, integrating equations (3.12) and (3.10) to obtain 

li = M(x a )e- T ~p , (3.37) 

and 

a 2 =Y?{x a )e- 2T (3.38) 

we can see that equations (3.6), (3.7) [(3.9) and (3.10)], and (1-12) [the analogue of (3.20)] are all 
invariant under 8^8, fi ^ p,, a ^ a where 8 2 = 8 2 + 'ip Ql p, — ^1 + p ,a 2 = a 2 . However, assuming 
that p > (i.e., positive pressure), we see from (3.36) that the model is only valid for r < r c , 
for some critical value r c (so that 8 2 > 0), and for early times (r ^ -00 or t ^ + ) the model is 
indistinguishable from the zero-pressure model. 
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B. Pressure-free case. 

The special case of dust is of particular importance since in this case the constraints from the 
imposition of self-similarity are the least restrictive and also because there are examples of dust 
models admitting a kinematic self-similarity (Carter and Henriksen, 1989). We also note that in the 
case of dust (p = 0), it follows that C^T a i, = (a + 2a)T a f } ; that is the total energy-momentum tensor 
satisfies an equation of the form (1.12). 

When p = 0, the differential equation (3.23) reduces to 



with the solution 
where (erf 7^ 0) 



ip" + 3V>' + 2-0 = 0, 
iP = 6 2 = c{x a )e- 2T + d{x a )e- T 
1/2 2V~c 



t 



¥H4 



(where the function of integration is chosen to ensure that t — > + as r 
can write 



4 t 



2^~c 



Ac 



2^1 



where 

Integrating equation (3.28) yields 
and equation (3.25) gives 



2j~c 



c 
d' 



fi = M(x a )e- T , 



a = L [x )e , 
whence equations (3.6) and (3.7) (or (3.27)) are satisfied when 

M=\d, 
3 C - 

In addition, since a = —2a (p, 7^ 0), equation (1.12) yields 

M, a £ a 



2a + 



M 



T, a £ a ~ 0C = 0, 



where 



Using equations (3.42), (3.43) and (3.46), equation (3.48) reduces to 



t + 



t + 



2j~c 



(3.39) 
(3.40) 

(3.41) 

-00). Alternatively, we 
(3.42) 

(3.43) 

(3.44) 
(3.45) 

(3.46) 
(3.47) 

(3.48) 
(3.49) 

(3.50) 



and hence c = (and a is necessarily non-zero) and 
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Equations (3.14), (1.11), (3.9) and (3.11), remain to be satisfied; these equations (through equation 
(3.1) and due to the definitions of the shear and expansion in terms of the metric functions, their 
time derivatives, and the inverse metric functions), constrain the metric functions H a p. We recall 
that the resulting spacetimes are 3-Ricci flat (i.e., 3 R a b = 0). 

Let us investigate the asymptotic behaviour of the solutions. As t — > oo (r — ► oo), we find that 

9 = y/de~ T ' 2 = 2t-\ (3.52) 
and from equations (3.46) and (3.47) we have that 



£ - \ and J -> 0. (3.53) 



Therefore, the models are asymptotic (at late times) to an exact zero-pressure and flat FRW model 
with a power-law scale function; we note that this Einstein de Sitter model admits a homothetic 
vector. 

On the other hand, as t — > + (r — > — oo), we find that 

6» = Vce- T = i, (3.54) 
and from equations (3.46) and (3.47) we have that 

^ — and ^ — |. (3.55) 

Hence the models are asymptotic (at early times) to an exact vacuum, 3-Ricci flat solution; in 
particular, this exact solution has 

,i = 0, 

9=-: a 2 = -9 2 = \r 2 , 
f 3 3' 

3 Rab = 0. 

From equation (3.1) and from the definitions of the shear and expansion, these equations yield the 
following equations for the metric functions H a p(tx 7 ): 

H aP H a[j , t = - t , (3.56) 

H^H^H a p it H lS ,t = ^, (3-57) 

and (from equation (3.9)) 

H a (iM — H' 6 H a f,tHps,t + —H a p t t = 0. (3.58) 

In the special case of spatial homogeneity we necessarily obtain the Kasner model. In particular, 
in general in the case that H a/ 3(tx 7 ) is diagonal, i.e., H a p = diag {hi(tx 7 ), h 2 (tx 1 ), ^(tx 1 )}, 
these equations can be integrated to yield 

/ii = A 1 (x^)t 2ai , h 2 = A 2 (x^)t 2a2 ,h 3 = A 3 (x^)t 2a3 , (3.59) 

where the constants in (3.59) obey 

ai + «2 + a 3 = ai 2 + a 2 2 + a 3 2 = 1, (3.60) 

and hence we obtain the Kasner model. We note that this vacuum Bianchi I exact solution also 
admits a homothetic vector. 

Finally, we note that the spherically symmetric dust solutions of Lynden-Bell and Lemos (1988) 
and Carter and Henriksen (1989), that are a particular case of the perfect fluid solutions of Bcnoit 
and Coley (1996), represent a special non-trivial solution of the above equations. 
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C. Case of Zero Expansion. 

We shall also consider the special case in which 8 — 0, although this case is not of interest from 
a cosmological point of view. Let us choose comoving coordinates so that (with u = 0) 

ds 2 = -(U 2 )dt 2 + H al3 {tx~<)dx a dx 13 , (3.61) 

where 

u a = ^5g. (3.62) 

Using these coordinates we have that 

ii a = (0, [InUU), (3.63) 

6 = implies that 

H^Hau = 0, (3.64) 
and <I> = implies that <3? ;t = 0. In addition, from equations (1.10) and (1.11) we obtain 

f = (£°(*),r(* 7 ))- (3-65) 
From equations (2.31) and (2.39) we have that 

M,t = = <T, t , (3.66) 

and from equation (3.41) we find that 

(ti+p)[lnU], a = -p, a . (3.67) 

The class of static solutions is contained within this special case presently under consideration. 
Immediately, equation (2.36) yields 

{a + 2)n- {b + 2)p = 0, (3.68) 
whence equations (2.33) and (2.35) yield 

8(1 - a)a 2 + (a + 2)/x + 3(6 + 2)p = 0. (3.69) 
These two equations are best dealt with separately in two different subcases. 

Subcase (i): Either b = — 2 or p = 0. From equation (3.68) a = —2, and hence from equation 
(3.69), the shear is zero, a 2 = 0. Neglecting the case p = (since in this case we obtain ii a = from 
eqn. (2.32) and hence the resulting spacetime is a special static FRW spacetime), we consequently 
obtain 

a = b= -2, (3.70) 

and hence 

= 0; (3.71) 

i.e., 

H aPlt = (3.72) 

[and consequently H a p can be diagonalizcd— H a p = diag {/ii(x 7 ), h 2 {x 1 ), ft, 3 (x 7 )}]. From equations 
(2.11) and (2.33) we then obtain 

3 R = 2^(x 7 ), (3.73) 

and 

u a , a = ^ + -ip). (3.74) 
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In addition, from equation (2.9), we have that 

H a b = 0, 

i.e., the magnetic part of the Weyl tensor vanishes (note clumsy notation here), and equation 
yields 

E a b —^Rab — g 3 Rh a b- 

From equation (1.10) we have that 
and equations (1.12) yield 

Finally, equation (3.67) gives 

(ji+p){lnU) >a Z a = -p >a ?*, 
and (the integrability conditions of (3.67) give) 

P,a^,/3 = (*,aP,l3- 

From this equation we obtain the solution 

Inp = F(t, In^i), 

whence on defining T by 



equation (3.79) becomes 



with the solution 



which can be written as 



~dt 



dF dF 

+ 7^ = 1. 



d(lnn) dT 



F=-(lnn + T) + f(lnn-T), 



p = fj,f(lnn-T), 



where / is an arbitrary function of a single variable. Equations (3.67), (3.77) and (3.80) then yield 



(lnU), T =-($-a) 



f + f 



1 + / 



and 



which has the solution 



InU 



(IriU), 



/ + /' 



1 + / 



(Zn/i), 



f(m) + f'(m) 
1 + f(m) 



dm, 



(3.86) 



where m = (Zn/i — T). 
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Now, 

u a , a = H a P(x 5 )[(lnU), a p + (lnU), a (lnU),0 - T s a0 {x 5 )(lnU)^ (3.87) 
whence, on denning 

T{m) = J i + /m rfm ' (3 - 88) 



equation (3.74) becomes 



[T']{H^(ln^) Mj - H a erZ p (lnn)„} 

+ \T" - (T'f]{H a P{lnri, a {lnri,p} 

{fi} = 0. (3.89) 



+ 



|(l + 3.f) 



This equation contrains the metric functions H a p (x 7 ) and ^{x 1 ) [and /(Zn/i — T(t)) (J 7 ) and £°(t)]. 
The H a p [and £ a (x 7 )] are further constrained by equations (3.73), (3.75) and (3.76) [and eqns. 
(3.78) and (1.11)], where /z is given in terms of the metric functions through the EFEs. 

We notice that equation (3.89) is of the form of the sum of (three) terms in which each term in 
square brackets depends on t [through eqns. (3.65), (3.83), (3.85) and (3.88)] and each term in curly 
brackets does not. A particular solution of equation (3.89) has f = constant . Indeed, if / is not 
constant then in general (for \i ^ 0) we have that 

T" - {T' f = ci(l + 3/); T' = c 2 (l + 3/) [c 2 ± 0], (3.90) 

which on integration yields e - ^ = e?2 — ^d\e Cim l c ' 2 (ci ^ 0) or e _JF = d^ — d\m (c\ = 0) and hence a 
particular functional form for /(m), whence for consistency equation (3.88) implies that d\ = (in 
either case), which contradicts the assertion that / is not constant. 

Therefore, in general / is constant and hence from (3.85) we have that 

p = (7 - 1)M; / = 7 - 1 (const.). (3.91) 
In particular, we note that p = p(x 7 ). Finally, equation (3.67) gives 

[lnU], a = (j-^j (lnti),a, (3.92) 

which integrates to yield 

U = u(t) f i- 1+1 ^. (3.93) 

However, since U is separable, a redefinition of the time coordinate in (3.61) can be employed to set 
u(t) = 1. Hence 

U = t/(x 7 ) = [i~ 1+lh , (3.94) 

and the metric is independent of time and hence the model is completely static. 
Equation (3.77) then yields 

= (a+ 2(1 ~ 7) ^)t + c, (3.95) 

whence the £"(x 7 ) are constrained by equations (3.78) and (1.11) while the i? a( g(x 7 ) are themselves 
subject to equations (3.73), (3.75) and (3.76). 
Finally, from (3.91) we have that 

F(m) = ( ^^-dm = ^—^-(ln^-T), (3.96) 
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and equation (3.86) yields 



lnU = ~ J(a-t t )dT-T{m) 



Q — — -(lnp-T) = - — -Infi + d, (3.97) 

7 7 7 

which is consistent with equation (3.93). In addition, in this case equation (3.89) reduces to 

7 

+ H^Tl (ln^ n = ^"^jV (3-98) 

This equation must be satisfied in addition to equations (3.73), (3.75) and (3.76). 
Subcase (ii): In this subcase we have that 6^—2 and p ^ 0, so that equation (3.68) yields 
p = wnence equations (1.12) imply that a = b(^= —2) and hence p = fi (note that p_ t = 0), 

and finally equation (3.69) yields 2(1 — a)a 2 + (a + 2)p = 0; that is, we have that 

(b + 2)p ± 0, a = b, p = p, o 2 = M- (3.99) 

From equations (2.11) and (2.33) we also obtain 

3 i?=-ff±^ = ^V (3.100) 
(1 - a) 

From equation (3.67) we obtain 

U = u(t)fi-^ (3.101) 

(which we note is the special case of (3.93) with 7 = 2— corresponding to p — p), whence by 
redefining the t-coordinates to set u(t) = 1 we have that 

U = pT*. (3.102) 

Therefore, although p,p and U are independent of t, the metric functions H a p = Ha^tx 1 ) depend 
on t and the models are not static (since a 2 ^ for [i 7^ 0, cr a h 7^ and hence H a p^ cannot vanish 
for all afl). However, the H a p are constrained by equations (3.64), (3.99) and (3.100), in addition 
to equations (2.34) and (2.37), etc. For example, using (3.63) and (3.64), equation (3.100) yields 



U[U aP - T s ap U. s ]H afi = l a + 2 ^ (const.). (3.103) 

(1 - a) 



Using equations (3.99) and (3.102), equations (1.12) yield 



(lnU), a C = -\= —{InnlaC, (3-104) 



whence equation (1.10) yields 
which integrates to give 
Equation (1.11) remains to be satisfied. 



£,°(t) = ^(a + 2a)t + c. (3.106) 
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4. More Special Cases 

There are two special cases of particular interest in which the kinematic self-similarity is either 
parallel or orthogonal to the velocity vector. Let us next consider these two cases separately. 

A. £ parallel to u. 

In Coley (1991) it was shown that if a perfect fluid spacetime admits a proper HV parallel to the 
velocity vector then that spacetime is necessarily an FRW spacetimes satisfying 9 = — ^9 2 ,9 ^ 
[with the very special and physically unreasonable equation of state /z + 3p = 0; strictly speaking 
this result was only proven in the case p = p(/i) and \i + p ^ 0] . 

Suppose that a spacetime admits a kinematic self-similar vector £ parallel to u, i.e., 

C = Au a , (4.1) 

then, from equation (2.12), 

A(u a . b + u b . a ) + A^ b u a + A ta u b = 2g ab + 2(1 - a)u a u b . (4.2) 
Contracting equation (4.2) with g ab and u a u h in turn then yields 

A = a, AO = 3, (4.3) 

and hence 

9 + ^9 2 =0. (4.4) 
Contracting equation (4.2) with h c a h b d then yields 

(Tab = 0, (4.5) 
i.e., the spacetime is shear-free, and finally contracting equation (4.2) with u a yields 

u a = (InA- 1 )^. (4.6) 
Assuming that the vorticity is zero, equation (2.6) then yields 

9,bh h c = 0. (4.7) 
Therefore, from equations (4.3) [A = 3/9] and (4.7) we see that 

u a = 0, (4.8) 

i.e., the acceleration is zero. 

Hence, since the shear, vorticity and acceleration are all zero, the spacetime is necessarily FRW 
(Ellis, 1971; Coley and McManus, 1994). Therefore, there exists coordinates in which the metric is 
of the form (4.48) and equation (4.4) yields 

9 t + ^9 2 ^0; 8^-r 1 (4.9) 
3 a 

(defining the constant of integration so that 6 — > oo as t — > + ). The remaining equations to be 
solved reduce to 

a + 2a 3(6 + 2a) 

rP = °> ( 4 - 10 ) 



2(1 -ay 2(1 -a) 
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3(a + 2) 3(6 + 2) 
2(1 - a)^ ~ 2(1 - a)' 
29 2 (l - a) = -3fi- 9p, (4.12) 



fi--) = 4^-T±% (4.1D 

V 3/ 2(1 -aY 2 1-af' v ' 
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and 

fi+(ji + p)0 = O. (4.13) 
It can be easily shown that a ^ —2a and p = lead to contradictions, so that (4.10) yields 

a = b = -2a. (4.14) 

Equations (4.11)-(4.13) then yield the consistent solution 

'2a 



P = 

and 



3 ^ 



ix (4.15) 



»=-/ = - 2 t-\ (4.16) 
3 or 

We note that the above solution reduces to the homothctic case when a — 1 (see, e.g., cqns. 
(4.4) and (4.15) and (4.16)). However, unlike the homothetic case there exist models with realistic 
equations of state; for example, < p < fi for | < a < 3. 

B. £ orthogonal to u. 

Mcintosh (1975) showed that a perfect fluid spacetime cannot admit a non-trivial homothetic 
vector which is orthogonal to the fluid 4-velocity unless p = [i. 

We shall work in comoving coordinates in which the metric is given by (3.61) and the velocity 
vector by (3.62). In these coordinates the acceleration is then given by 

u = 0;u a = (lnU) >a . (4-17) 

Since £ is orthogonal to u, £° = 0, and equation (1.11) implies that £ Q = £ Q (x 7 ), and consequently 
equation (1.10) reduces to 

U^ a = aU. (4.18) 

Equations (1.12) yield 

/V£ Q = W (4-19) 

and 

P,*Z a - bp. (4.20) 
Also, we can further specify the coordinates so that 

?*=Z{xi)52, (4.21) 

and in these coordinates (if aab ^ 0) equations (4.18) - (4.20) yield 

(in[tf 1/a D,* = {lnW a ]) /x = (MP 17 "]),* = |, (4-22) 
which can be partially integrated to yield 

H - 7(*. V, z)U a/a ; P - g(t, y, z)U h ' a . (4.23) 
Finally, the conservation equation (2.32) reduces to 

(p + p)(lnU), a = -p, a , (4.24) 
whence, on contraction with £ a and using equations (4.18) and (4.20), we obtain 

a\n + p] = -bp. (4.25) 
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Apart from the special subcase a = b = 0, which implies that 

U x = 0=p, x , (4.26) 

equation (4.25) implies that (i) if b = 0, then either a = (the special subcase above) or \i + p = 
(a case which we shall not consider here), (ii) if a = 0, then either b = (special subcase) or p = 
(a case considered earlier), and (hi) if b = —a, then either a = (and b = 0) or ji = (a case of no 
interest here; however, see the appendix), otherwise we have that in the general case 

p=—^-^ ab(a + b)^0. (4.27) 

(a + b) 

Therefore, equations (4.27) and (1.12) immediately imply that 

a = b. (4.28) 

The conservation equation (2.31) then yields 

fj, tt = -a(a + a)~ 1 fj,Ue, (4.29) 

which can be regarded as an equation for 6 in terms of \i and U. 

Using equations (4.27) and (4.28), equations (2.33) and (2.35) then yield 

e + ^ = -2o>- (; f +2 X°- 2a K , (4.30) 



equation (2.36) yields 



a2 _ 3(a + 2)(a + 2a) 



from which we deduce that 



9 + 9 = T7^ T? , > » ( 4 - 31 ) 

4(1 - a) (a + a) 

\e 2 = 2a 2 + £±*L. (4.32) 
3 (1 — a) 

Differentiating this equation with respect to t, and using equations (4.29) and (4.31), we obtain 
(a ? 0) 

a = -<t(9, (4.33) 
which also follows from equation (2.37). Finally, from equations (4.32) and (2.45) we obtain 

t„ (a + 2a) 

* = -W^ (4 - 34) 

Now, defining F(x,y,z) by {lnF) tX = l/£, equations (4.22) then yield 

/i = /(t,y,z)F a , (4.35) 

U = g(t,y,z)F a . (4.36) 



Since F it = 0, equation (4.29) yields 



a fg 



whence equation (4.31) then yields 



fg 



(a + a) (f, t ) 2 3a(o + 2)(o + 2a) +2a 
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Because the only term in (4.38) that depends on x is the term F a+2a (and since F tX ^ from 
(lnF) iX = l/£ 7^ 0), and assuming that a ^ (since a = implies 6 = 0), the only way that this 
equation can then be satisfied is for either 

(i) a = -2, or (ii) a = -2a. (4.39) 

Hence, from equation (4.31) we have that 

9 + 9 2 = 0. (4.40) 
From equations (4.38) and (4.40) we can then integrate to obtain 

g = J{ yi z)r {a+2a)la f,t. (4.41) 

Finally, let us consider the two cases in (4.39) separately. 

(i) a = —2. In this case (4.25) implies that 

P=^ (4.42) 

and equations (4.32) and (4.34) imply that 

a 2 = V, (4.43) 

and 

3 R = 2/i. (4.44) 

(ii) a = —2a. In this case (4.25) implies that 

p = n (4.45) 
(i.e., the fluid is stiff), and from equations (4.32) and (4.34) we have that 

H + cj 2 = 1 -8 2 , (4.46) 

and hence 

3 R = 0. (4.47) 

Further progress (e.g., further constraining the form of the functions F(x,y, z), f(t,y, z) and 
J(y,z) or determining the form of the metric functions H a p(t,x,y, z)) can be made in specific 
spacetimes. For example, spherically symmetric spacetimes have been studied [and it has been 
claimed that in such spacetimes if there exists a vector field £ satisfying (1.11) which is orthogonal 
to u, then the resulting spacetime metric is singular (Ponce de Leon, 1993)]. We shall investigate 
the existence of kinematic self-similar vectors in the special case of FRW spacetimes in the next 
subsection. 

C. FRW models. 

We recall from Maartens and Maharaj (1986) that (perfect fluid) FRW models, when written in 
the comoving form 

dr 2 



ds 2 = -dt z + R A (t) | - _ ^ 2 + rW j , (4.48) 

admit a homothctic vector P G (= dt) parallel to u for all k when R(t) = dt, d constant (and hence 
/i = —3p = 3(1 + kd~ 2 )t~ 2 -. [i + 3p = 0), and admit a homothetic vector H (= td t + rd r ) when 
R = dt c but only in the case k = (whence /i = 3c 2 i~ 2 , p = (7 — l)fi with 7 = 2/3c). 
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We cannot simply apply the results from section 3 here since it was explicitly assumed there that 
(7^0. However, when a = lo = u a — 0, we immediately obtain from equation (2.35) that 

(a + 2a)n + 3(b+2a)p = 0. (4.49) 

This again suggests two subcases. First, when a = —2a, we have that either b = —2a or p — 
0, whence from equations (2.33) and (2.36) we obtain (in either case) /j, = ^6 2 and hence from 
equation (2.11) we obtain 3 R = 0. Second, if a ^ —2a, then if b = a, from (4.49) we obtain 
/i + 3p = [on the other hand, if we assume a ^ d ^ —2a (p ^ 0), then from (4.49) we have that 
p = [—{a + 2a)/'3(b + 2a)]/i, whence from equations (1.12) we deduce that a = b, resulting in a 
contradiction]. Therefore, in order for an FRW spacetime to admit a kinematic self-similarity then 
necessarily either k = (zero-curvature) or fi + 3p = 0, in direct analogy with the homothetic case. 

However, we shall not proceed with this type of analysis here since we wish to determine and 
display all the kinematical self-similarities in conventional forms; i.e., in the form corresponding to 
(4.48) and in a form that explicitly displays the self-similar form of the solutions. 

In order to derive the spacetimes in manifestly self-similar form, we write the metric in comoving 
spherically symmetric coordinates adapted to £, viz., 

ds 2 = _ e 2m dt * + e 2<M€) dr 2 + r 2 5 2(^ dn 2 i (4 5Q) 

where dfl 2 = d6\ + sin 2 d\dB\, and 

» = e-*£, (4-51) 

where, assuming £ is neither parallel to nor orthogonal to u, the kinematic self-similar vector £ and 

the corresponding self-similar variable £ can be written in one of the following forms (Carter and 
Henriksen, 1989): 

first kind (a = 1; HV) : £ = t— + r— ; £ = r/t, (4.52) 

at or 

second kind (a ^ 0, 1) : £,=at—+r—; i = r/(at) x / a , (4.53) 

at or 

d d 

zeroth kind (a = 0) : £ = — + r— ; £ = re"* . (4.54) 
v ' dt dr 

Now, if the acceleration is zero, we have that = ao, a constant, and if the shear is zero, we 
have that = f Q S, where f a is a constant. Therefore, the FRW metric can be written in the form 

ds 2 = -a 2 D dt 2 + S 2 (£)[f 2 dr 2 + r 2 dn 2 }. (4.55) 

(i) If a ^ 0, and assuming S' ^ (non-static case), the EFEs (for a perfect fluid source) then yield 

S = d£ c , (4.56) 

and (c^ —1) 

f 2 = (l + c) 2 . (4.57) 

The forms for £ in (4.52) and (4.53) are invariant under the changes t — > at and r — > br, and a 
and b can be chosen to set a a = 1 and d — 1, and we can then write the metric as 

ds 2 = -dt 2 + e c {j^-y +rWj. (4.58) 

Recalling that £, — r{at)" 1 / a {a ^ 0), and defining a new radial coordinate f by f = a a r c+1 , where 
a = —c/a, the FRW metric takes on its familiar form 



ds 2 = -df 2 + t 2a {df 2 + f 2 dfl 2 }, 



(4.59) 
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where £ is now given by 

d d 

t = ut- + {l + c)r-_. (4.60) 

We note that each orthogonal 3-space is flat (k — 0) and that the simple power-law solutions 
give rise to an equation of state of the form p = (7 — and equations (1.12) are automatically 
satisfied. The known homothetic case (£ = H for k = 0) is included here. In fact, each spacetimc 
(4.59) (for each value of a) admits kincmatical self-similar vectors for all values of a (7^ 0). [This is 
known to be true in the dust subcase (Carter and Henrikscn, 1989)]. This, of course, arises due to 
the separability of 5(£) in (4.55) (e.g., see equations (4.53) and (4.56)). 
(ii) If a = 0, then the EFEs imply that 

S = cf, (4.61) 

and 

.fo 2 =4. (4.62) 
Since £ = re~*, under the transformation t — > t + b, b can be chosen so that c = 1, and hence 

ds 2 = -a 2 a dt 2 + r 2 e- 2t (4dr 2 + rW). (4.63) 

Finally, defining the new variables i = a a t, f = r 2 , the metric becomes 

ds 2 = -dP + e- 2t ' a °{df 2 + f 2 dn 2 ), (4.64) 

and the kinematical self-similar vector is given by 

* = a °m +2f w (4 - 65) 

Again, this FRW model is flat (k — 0). However, this case has no homothetic vector analogue. The 
equation of state for this FRW model is given by \i = — p = I2a~ 2 , and the metric is, of course, the 
de Sitter metric. Hence, the de Sitter model admits a self-similarity of the zeroth kind. Equations 
(1.12) arc trivially satisfied with a = b = 0. 

Finally, we must deal with the two special cases not considered above in which £ is either parallel 
to or orthogonal to u. For illustrative purposes let us return to the more conventional coordinate 
system in which the metric is given by (4.48). Writing £ in the form 

^ = rM) l +rM) ^' (4 - 66) 



equations (1.10) and (1.11) yield 



where c is an arbitrary constant, and 



and 



C - at + /3, (4.67) 
f = ( 2— Sir, (4.68) 



2 R 

— (at + /3) = c, (4.69) 



fc(c-2) = 0. (4.70) 



From this last equation we see that either k = 0, and the FRW spacetime is flat, or c = 2, whence 
the kinematic self-similar vector is parallel to the fluid velocity vector (i.e., £ r = 0). Let us consider 
this latter case first. 
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(a) £ parallel to u. In this case c = 2 (i.e., £ r = 0) and k is unrestricted by equations (4.67)-(4.70), 
whence equation (4.69) becomes (a 2 + /3 2 ^ 0) 

|(at + /?) = l. (4.71) 
If a ^ 0, we can set /3 = by a time translation, whence equation (4.71) yields 

R = dt 1 ' a 1 (4.72) 

so that when a = 1 (homothetic vector case) we recover the usual FRW model with \i + 3p = 
(valid for all fc), whereas for a ^ 1 all power law solutions of the form (4.72) admit a vector field 
satisfying equations (1.10) and (1.11). 

We note that if k = 0, then the spacetime (4.59) also admits additional kinematic self-similarities 
parallel to u (in addition to those given by (4.60)). If, on the other hand, k ^ 0, we can see that 
the resulting FRW models with power law scale factors always admit a vector field parallel to the 
fluid velocity vector that satisfies equations (1.10) and (1.11). However, equations (1.12) need not 
be satisfied. Indeed, from equation (4.72) we find that 

M=^" 2 + ^- 2/a , (4-73) 
P=-(2-3/a)t- 2 -^t- 2 / a , (4.74) 

so that equations (1.12) can only be satisfied if (cither) a — 1, whence £ is a homothetic vector and 
M + 3p = as before (or k = 0). 

It is curious to note, however, that if we consider the perfect fluid source to be due to two separate 
comoving perfect fluids (Coley and Tupper, 1986), so that 



M = Mi + M2, 

P = Pl +P2, 



(4.75) 



where 



3 



a 

and 



Mi = ^t- 2 - Pl = (2a/3 - 1)mi (4.76) 



M2 = ^t- 2 /«; p2 = -\l>2 (4-77) 



then 



£«n^i = -2a^i, £ (u p 1 = -2ap 1 

£ ?|| M2 = -2M2, C i]]P2 = -2p 2 (4.78) 

where = a t§[, so that each separate fluid satisfies equations of the form (1.12). 
If a = 0, then equation (4.71) yields 

R = de 1 ! . (4.79) 

That is, FRW spacetimes with (4.79) admit a vector field £|| = cd§^ = atu satisfying equations 
(1.10) and (1.11). If k = 0, then spacetimes of the form (4.64) admit additional kinematic self- 
similarities parallel to u. There is no analogous result in the case of a homothetic vector. If k ^ 0, 
however, then equations (1.12) cannot be satisfied. Again, if we consider the two separate comoving 
perfect fluids interpretation (4.75) with 



_3_ 

J 2 
3k 



Mi = = "Pi' ( 48 °) 



M2 = ^e- 2t//3 = -Spa, (4.81) 
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then each separate fluid satisfies equations of the form (1.12). 

If c 7^ 2 (£ r 0), then from equation (4.70) we must have k — 0; i.e., the FRW models are of 
zero curvature. If a 2 + f3 2 ^ 0, then solving equation (4.69) yields the spacetimes (4.59) and (4.64) 
obtained earlier. If a = j3 = 0, then £° = 0. 

(b) £ orthogonal to u. If a = (3 = 0, then £° = and the vector £ is orthogonal to the fluid velocity 
vector. Equation (4.69) simply yields c — 0, so that 

that is, every flat FRW model admits a vector field of the form (4.82) which satisfies the conditions 
(1.10) and (1.11); in fact, 

C u u a = 0. (4.83) 
In addition, since fi — p(t) and p = p(t), it follows immediately that 

£ U fi = = £ u p, (4.84) 

so that equations (1-12) are trivially satisfied. Hence, every flat FRW model (for any R(t)) admits 
a vector field given by (4.82), orthogonal to u, that satisfies equations (1.10), (1.11), (1-12); in 
particular, £^ ± h a b = 2h a b and equations (4.83) and (4.84) are satisfied. 

5. Conclusion 

After a brief review of self-similarity and its applications in general relativity, the covariant notion 
of a kinematic self-similarity in the context of relativistic fluid mechanics was introduced (Carter 
and Hcnrikscn, 1989), which it was argued is the natural relativistic counterpart of self-similarity 
of the more general second (or zeroth) kind and hence is a generalization of a homothety which 
corresponds to self-similarity of the first kind (Cahill and Taub, 1971). Various mathematical and 
physical properties of spacetimes admitting a kinematic self-similarity were discussed. We note 
the relationship between a kinematic self-similarity and what Collins and Szafron (1979) term an 
intrinsic symmetry and what Tomita (1981) refers to as a partial homothety (however, see also 
Tomita and Jantzen, 1983). 

The governing equations (adopted from Ellis, 1971) of the perfect fluid (cosmological) models 
under investigation were introduced, and a set of integrability conditions for the existence of a 
proper kinematic self-similarity in the spacetime models was derived. These important constraints, 
the integrability conditions, given by equations (2.16) in general and by equations (2.27) - (2.30) 
in the particular case of a perfect fluid, played a central role in the resulting analysis. All of the 
relevant equations were then given in the physically important case of zero vorticity. 
A. Summary 

Exact solutions of the irrotational perfect fluid Einstein field equations admitting a kinematic 
self-similarity were then sought in a number of special cases. Since the integrability conditions 
constitute very severe constraints, such solutions are necessarily of a particularly simple form. 

First the geodesic (i.e., zero acceleration) case was considered in subsection 3. A. It was proven that 
(provided the shear is non-zero) the 3-spaces orthogonal to u are Ricci-flat, that is, 3 R a b = 0. This 
case consequently merits further consideration since there have been various studies of spacetimes 
with vanishing 3-Ricci tensor (see, for example, Collins and Szafron, 1979, and Stephani and Wolf, 
1986). It was further proven that a = —2a, and the form of the kinematic self-similar vector £ 
and the pressure were given in the particular coordinates (3.1)/(3.2) by equations (3.16) and (3.18), 
respectively (a = b = and p = constant in the special case a = 0). The expansion was shown to be 
governed by the differential equation (3.23)/(3.24), which was in fact integrated in the special case 
a = [see equation (3.36)]. 

The further specialization to dust (i.e., zero pressure) was then considered in subsection 3.B. In 
this case the governing differential equation (3.23) reduces to equation (3.32) which can be completely 
integrated to obtain the expansion (see eqns. (3.40) - (3.43)). The asymptotic properties of this 
class of solutions was studied, and it was found that the resulting models are asymptotic (at late 
times) to an exact flat FRW (Einstein-de Sitter) model and are asymptotic (at early times) to an 
exact vacuum (3-Ricci flat) model with 3er 2 = 9 2 = t~ 2 — for a large class of models this was shown 
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to be the Kasncr (Bianchi I) solution. We note that these exact (asymptotic) solutions are known 
to admit a homothetic vector. 

The case of zero expansion, studied in subsection 3.C, was shown to subdivide into two subcases. 
In the first subcase the models are necessarily shear-free with zero magnetic part of the Weyl tensor, 
and in general p = (7 — l)fi (and equation (3.94) is satisfied) and the resulting models are completely 
static. In the second subcase the models are necessarily stiff (p = ji) non-static perfect fluid models. 
In the coordinates (3.61)/(3.62) the form of the kinematic self-similarity was found to be given by 
equation (3.65) and either equation (3.95) or (3.106). 

In subsection 4. A perfect fluid spacetimes admitting a kinematic self-similar vector £ parallel to 
the velocity vector u were studied, and it was proven that such spacetimes are necessarily FRW 
spacetimes with p = — and \i = \Q 2 = -§it~ 2 . In the case that £ is orthogonal to u, it was 
shown in subsection 4.B that in general [ab(a + b) 7^ 0] p = (7 — l)/i with 7 = b/(a + b) and 9 
satisfies the differential equation (4.40) (and a = b), and that either a = — 2 whence equations (4.42) 
- (4.44) follow or a — -2a and equations (4.45) - (4.47) result. 

Finally, in subsection 4.C, the existence of kinematic self-similarities were studied in FRW space- 
times. In the general case in which £ is neither parallel nor orthogonal to u, it was shown that if 
a =/= then the resulting FRW model is flat with p = (7 — l)fi and the scale factor is of a simple 
power-law form [and £ is given by (4.60) in conventional coordinates]. Indeed, each such FRW model 
admits kinematic self-similar vectors for all values of a 0) in addition to a homothetic vector. On 
the other hand, if a = the FRW model is necessarily a flat de Sitter model with \i + p = [and £ 
is given by (4.65)]. Hence the de Sitter spacetime admits a self-similarity of the zeroth kind (note 
that such a spacetime cannot admit a homothety). In addition, all FRW models with a scale factor 
of the power-law form (4.72) with a / [and with a scale factor of the exponential form (4.79) 
when a = 0] were shown to admit a vector field £ parallel to u which satisfies equations (1.10) and 
(1.11). Here the curvature need not be zero; however, in the flat case the FRW models admit these 
special vector fields (parallel to u) in addition to the kinematic self-similar vectors mentioned above. 
Finally, every flat FRW model (with a scale factor of any form) was shown to admit a vector field £ 
of the form (4.82), orthogonal to u, which satisfies equation (1.10) and equations (1.11) and (1.12) 
[see equations (4.83) and (4.84)]. 
B. Discussion 

There are a variety of circumstances in general relativity theory, and particularly in cosmology, 
in which self-similar models act as asymptotic states of more general models. Indeed, in a number 
of classes of perfect fluid cosmological models with equation of state p = (7 — 1)^ and in which the 
governing equations reduce to a dynamical system, including, for example, spatially homogeneous 
models and silent universe models, and in some cases spherically symmetric models and G2 models, 
it is known that exact solutions admitting a homothetic vector play an important role in describing 
the asymptotic properties of these models (see Coley, 1996, for a review and appropriate references). 

For example, orthogonal spatially homogeneous models have attracted much attention since the 
governing equations reduce to a relativity simple finite dimensional system of autonomous ordinary 
differential equations (Wainwright and Ellis, 1996 - henceforward WE). Wainwright and collabora- 
tors (see Refs. in WE and Coley, 1996) have utilized an orthonormal frame approach and introduced 
an expansion-normalized (and hence dimensionless) set of variables to study these models. In partic- 
ular, it was proven that all the singular points of the (corresponding "reduced" ) system (of ordinary 
differential equations) correspond to exact solutions admitting a homothetic vector (Hsu and Wain- 
wright, 1986). It is in this sense that self-similar models play an important role in describing the 
dynamics of spatially homogeneous models asymptotically. [The situation is complicated by the fact 
that in the more general classes of Bianchi models there exist more complicated attractors than 
simple singular points; for example, in models of Bianchi type IX (and VIII) there is oscillatory 
behaviour with chaotic-like characteristics as one follows the evolution into the past towards the 
initial singularity due to the existence of a 2-dimensional attractor in the 5-dimensional phase space 
(WE).] 

In addition, in the class of inhomogeneous G2 cosmological models (in which the spacetime admits 
two commuting spacelike Killing vectors acting orthogonally transitively) it has been shown by 
Hewitt and Wainwright (1990) that the governing Einstein field equations can be written as an 
(infinite dimensional) autonomous system of first-order quasi-linear partial differential equations in 
terms of two independent dimensionless variables, and it was proven that the associated dynamical 
equilibrium states correspond to exact cosmological solutions that admit a homothetic vector (and 
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which arc consequently self-similar). In a particular subclass of G2 cosmologies, the separable 
diagonal G 2 models, it was shown that the models do indeed asymptote towards the dynamical 
equilibrium points, and Wainwright and Hewitt have conjectured that this may be the case for 
more general G2 models. Hence self-similar models may play an important role in describing the 
asymptotic dynamical behaviour of these inhomogeneous cosmological models. 

In this paper we have studied models which admit a kinematic self-similar vector. We note that 
in all cases in which we have either been able to integrate the equations to obtain exact solutions 
or we have been able to determine the asymptotic behaviour of a class of models, the asymptotic 
behaviour has been represented by an exact solution admitting a proper homothetic vector (and 
hence a model which is self-similar of the first kind). For example, in the pressure- free case studied 
in subsection 3.B, models were found to be asymptotic to the Einstein-de Sitter model (to the future) 
and the Kasner model (to the past), and both of these exact solutions admit a homothetic vector. In 
addition, the particular, exact, perfect fluid spherically symmetric solutions studied by Benoit and 
Coley (1996) [which include the dust solutions of Lyndcn-Bcll and Lemos (1988) and Carter and 
Henriksen (1989)] were shown to be asymptotic (both to the past and to the future) to exact FRW 
models which admit a homothetic vector, and it was also shown in Benoit and Coley (1996) that 
the equilibrium points at finite values of the autonomous system of ordinary differential equations 
which govern the general class of perfect fluid spherically symmetric, kinematic self-similar models 
correspond to exact solutions which admit a homothetic vector. It appears that the same is true for 
the analogous solutions in the case of plane symmetry (work in progress). 

It would be interesting to determine whether this is true in more generality. That is, it is an 
interesting and important question to determine the conditions under which models admitting a 
proper kinematic self- similarity are asymptotic to an exact homothetic solution. This might then 
shed light on when self-similar models of the second kind play an important role in determining the 
"intermediate asymptotic" behaviour of solutions (Barcnblatt and Zcldovich, 1972), and the role of 
generalized self-similar models in describing the asymptotic properties of models. 

The exception to the above is the de Sitter model. This flat FRW model with equation of state 
(i = —p = constant (or equivalently with a cosmological constant) does not admit a homothetic 
vector and is not asymptotic to a solution that does. Hence it cannot be true that kinematic 
self-similar models are asymptotic to exact homothetic solutions under all circumstances. 

The de Sitter model does, however, admit a self-similarity of the zeroth kind (see subsection 
4.C). It is curious to note that the de Sitter model acts as an asymptotic state of more general 
models (according to the so-called cosmic no- hair theorems). In particular, it was shown by Wald 
(1983) that spatially homogeneous Bianchi models with a cosmological constant (except those of 
Bianchi type IX which recollapse) are future asymptotic to the de Sitter model. The existence of 
the de Sitter model as a future asymptotic state is indicative of exponential inflation. Hence, it 
is of interest to determine under what conditions the asymptotic states of cosmological models are 
represented by solutions of Einstein's field equations admitting a generalized self-similarity (i.e., not 
just a homothety). This question deserves to be studied further. 

Recently Wainwright 's work has been generalized to the case of imperfect fluid Bianchi models 
satisfying the non-causal linear Eckart theory of irreversible thermodynamics (Coley and van den 
Hoogcn, 1994) and the causal theory (both the truncated version and the full theory) of Israel and 
Stewart (see Coley et al., 1996, and references within). Dimensionless physical variables (similar 
to those used by Wainwright) were utilized and a set of "dimensionless equations of state" were 
assumed, whence it was again shown that in general the singular points of the resulting dynamical 
system are represented by exact homothetic solutions (Coley and van den Hoogen, 1994). In the 
exceptional cases the singular points correspond to models which violate the strong energy conditions 
and have constant expansion, and the models are analogues of the de Sitter solution (with the viscous 
terms mimicking a cosmological constant) and are consequently self-similar of the zeroth kind. 

Moreover, in other work viscous fluid models have been studied (particularly in the case of simple 
FRW and Bianchi spacetimes) in which the governing equations reduce to a (simple) system of 
autonomous ordinary differential equations, but since particular equations of state were assumed 
that are not of a "dimensionless" form the associated singular points do not necessarily correspond 
to exact solutions admitting a homothetic vector. [In this work the viscosity coefficients are modeled 
by both the non-causal Eckart and causal Israel-Stewart theories of irreversible thermodynamics; 
the reader is directed to the research papers of the Polish and Russian groups and the Russian and 
Spanish groups, respectively, which are fully referenced in the papers by Coley and collaborators 
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cited here.] In a preliminary investigation of this work it appears that all singular points correspond 
either to an exact solution which is known to admit a homothctic vector, to an exact zero-curvature 
FRW model (not necessarily admitting a homothety), or to a de Sitter- like solution. Hence, all 
of these models appear to be asymptotic to models which admit a kinematic self-similarity of the 
zeroth, first or second kind. Clearly this needs to be studied further. 

However, it is clear that kinematic self-similar models play an important role in describing the 
asymptotic properties of cosmological models. It is interesting to ask to what extent are cosmological 
models (or, rather, what is the class of solutions of Einstein's field equations which are) asymptotic 
to self-similar solutions, when self-similarity is understood in its more general sense. 

Appendix: The vacuum case. 

Vacuum spacetimes admitting a homothetic vector were studied by Mcintosh (1975), in which it 
was shown that a non-flat vacuum spacetime can only admit a non-trivial homothetic vector if that 
homothety is non-null and is not hypersurface orthogonal. 

The study of kinematic self-similar vectors in vacuum spacetimes is not physically well-motivated, 
since in general there does not exist a physically or intrinsically defined timelike vector u with respect 
to which the metric can be uniquely decomposed (unlike in the case of a perfect fluid spacetime, for 
example, in which there exists such a vector which has both an intrinsic physical and geometrical role 
- it is both tangent to the fluid flow and is the unique normalized timelike eigenvector of the Ricci 
tensor), and consequently there exists no such u with respect to which the definition (1.10) /(l.ll) 
can be applied. However, for curiosities sake, let us study the consequences of the existence of 
an intrinsically defined timelike vector u which satisfies equations (1.10) and (1.11) in a vacuum 
spacetime. 

From subsection 2B (for ui = 0), in the case /i = p = we obtain the equations 

W C ; C = 0, (A.l) 

and 

6 + 6 2 = 0, (A.2) 

and hence 

a 2 = V, (A.3) 

and consequently 

3 R = 0. (A.4) 

Let us adopt coordinates so that the metric functions are defined through (3.61) and equations 
(3.62) and (3.63) are valid, whence on defining the new time parameter, r, by 

^ft^OTtt' (A ' 5) 

we can then integrate equation (A.2) to obtain 

6 = e(x-')e- T , (A.6) 

whence equation (A.3) yields 

CT 2 = Ie 2 e- 2r , (A.7) 

o 

which also follows from (2.37) [when a ^ 0; if a = 0, then if follows that 9 = 0]. By definition we 
have that 

H afj H' af} = 2, (A.8) 



and 

2 3 



O a p = ^ H Lp - ^ H af3, (A.9) 
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and 

3 RaP = ^U, a p-^TZ fi U t ^ (A.10) 

(all other components of o a b and 3 R a b vanish). Finally, equation (2.37) yields 

a' a0 - H^a l( Ji' !i)a + a a0 = 0, (A.11) 

and using (A. 9) we consequently obtain 

H'^-H^H' ai H' 0S = O. (A. 12) 

In the case that H a p is diagonal, i.e., H a p(r 1 x' 1 ) = diag{/ii, hi, hs}, we can integrate equations 
(A. 12) to obtain 

K(T,xi) = G u ( X i)e 2F »^ T (i/ = 1,2,3), (A.13) 
whence equation (A. 8) then implies that 

F 1 +F 2 +F 3 = l. (A.14) 

We note that the spatially homogeneous vacuum Bianchi I (Kasner) model is a particular solution 
of equations (A.2) - (A.4), (A.8) and (A.12) in which U = 1 (and hence 3 R ab = 0). 
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